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ABSTRACT 

The squeezed states or Bogoliubov transformations and wavelets are applied to two prob- 
lems in non-relativistic statistical mechanics: the dielectric response of liquid water, e(q, u>), 
and the bubble formation in water during insonnification. The wavelets are special phase- 
space windows which cover the domain and range of L'nL 2 of classical causal, finite energy 
solutions. The multi-resolution of discrete wavelets in phase space gives a decompostion 
into regions of time and scales of frequency thereby allowing the renormalization group 
to be applied to new systems in addition to the tired ‘usual suspects’ of the lsing models 
and lattice gasses. The Bogoliubov transformation: squeeze transformation is applied to 
the dipolaron collective mode in water and to the gas produced by the explosive cavitation 
process in bubble formation. 

1. INTRODUCTION 

Water is extremely important in chemistry and biology both as a solvent and as a neu- 
tral medium which is rich in possibilities. 1 - 19 The dielectric response of a medium, e(q,u) 
is related to the index of refraction as, n(q,u>) = y/fjl, whose real part controls wave prop- 
agation and whose imaginary part gives the attenuation. An August 1993 dissertation at 
Missouri University by one of us, S.-H. Kim 19 , formulates a new model which is rather 
successful in accounting for the complex index of refraction n(u) + ia(u>), over 14 to 15 
decades of frequency. This model is a number- conserving relaxation-time, 7 Boltzmann 
transport equation model in phase space 5,6,9 based on Kubo’s fluctuation-dissipation the- 
orem 12 and a Kerr approximation. 11 The free rotation gives an IR peak at the correct 
frequency which is much larger than the one found experimentally by Simpson et al. 15 A 
phase-space analysis of our model based on wavelets will be presented here. 

The problem of bubble formation in liquids is both old and difiicult. 20-25 It is now 
known that two different mechanisms of bubble formation exist: cavitation, where the 
growth is explosive and nucleation where diffusion of molecules yields a slow growth. 
This work will be restricted to cavitation. All bubbles are interesting for electromagnetic 
wave propagation in water because they give additional spatial dispersion ^-dependence to 
the index of refraction of the medium through their Mie scattering. 

The approach taken here will be to form a mean-field theory of a theory in 4 space- 
time dimensions. It is a variant of an idea due to Kaup, 25 Wilhelmsson 26 and Glimm 27 
that bubbles are non-linear, coherent effects in fluids. Kaup’s clever title is, “Cavitons 
are solitons” for a plasma driven by both sound and electromagnetism. Our equations are 
completely different and have solitary wave solutions rather than solitons. These solitary 
waves will have finite lifetimes in contrast to the infinite lifetimes of Kaup’s cavitons, which 
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he first showed were an integrable system. However, since both models are subject to other 
perturbations, this difference is not actually as serious as it first appears. 29 ^ 

There is also a practical reason for trying to understand bubble formation. Apfel ' , 
Nybor g 32 , Rose and Goldberg 33 have shown that ultrasound at the frequencies used for 
medical diagnostics forms bubbles in water and tissue. They have raised questions about 
the safety of ultrasound but have not yet shown any evidence that these bubbles (or any 
other effect of ultrasound) are harmful. Trevena 34 has suggested that bubble collapse may 

occur through shock waves, which could damage cells. 

Both of these problems will be analyzed using wavelets 3 ®' 51 which are a special class of 

windows 

(i) which cover the domain ft 1 or 3ft 2 in refs. (36-49) and the Lorentzian manifold of 
special relativity in ref (50,51); and 

(ii) are dense in the range of L 2, D L x . 

The first space L 2 is required classically for finite energy signals and for a probabili- 
ty interpretation of Hilbert space states in non-relativistic quantum mechanics, and the 

space L 1 is required classically for causality and is a major technical requirement in non- 

relativistic quantum mechanics which provides absolute continuity. Clearly, (i) assures the 
existence of measurability of the underlying space and (ii) places additional restrictions on 
the solutions for the correct dynamics and symmetries. In particular, the symmetry group 
is the set of dilations. In the continuous one- dimensional case dilations are generated by 
the transformation D 

(D a b)(x) = ax + b C 1 ) 

(. D ab xp){x ) = |a| 1 / 2 V> (a® “ t>) ( 2 ) 

on the coordinates x and solutions xp 6 L 2 H L 1 , where a € 3ft+ and b € 3ft 1 . In the discrete 
case (k,j) € (the integers) with a = 2 (one choice, but not the only one) 

(D jk )(x) = 2 i-x + k ( 3 ) 

(Dj.*X*) = • >M*) = |2p* (** - *) • W 


In the discrete case, I 2 -> t 2 and L x - t x and xp kj “live in” i 2 D i x . The choice of a = 2 
gives the Lebeque covering of 3ft 1 and extends to 3ft" and thus any other integer with this 
property will suffice. To satisfy (2) the conditions (2a) completeness 


Vxp jk = L 2 C\L X 

(5) 

(2b) mean zero for 

ip, mean one for (p 



xp( 0) = f rp(x)dx = 0 
J— oo 

(6) 


- f°° 

m i p(x)dx = 1 

J — OO 

(7) 
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where the multiresolution decompostion of Mallet satisfies Vj = V hl WWj-u i.e. Vj.i U 
= Yi Vj-i ^ Wj. 1 = {0}. For fixed j , eqs. (3,4) reduce to a translation on t? . 

t 2 nt' 



If k is considered as a time (space) translation, kto, and V as a frequency (wavenumber) 
scale, 2 } u 0 , with t^up > 1 then each fixed jo corresponds to a discrete translation group 
at frequency scale 2^°u>o and changing jo changes the frequency scale. Even if tj) or <f> axe 
time (space) signals, the dilation operators D ab or D jk map them into the time (space)— 
frequency (wave number) phase space. Thus 

<t>{x) € Vj and i>(x) e Wj (8) 

implies that 

<f>(2x) € Vj . i and ip(2x) e Wj . x . (9) 

The interpretation of the V"s spanned by the <f>y s is as a low band-pass window at frequency 
scale 2 j u>o whereas the higher frequencies are in the Wj s spanned by ipjjs at that frequency 
scale. 

From this short discussion two things should be obvious. One is that although a wavelet 
is a window, it is a great deal more as it must cover each space-time point in order 
not to miss any point particles, or sources of charge or mass. The other is that all of the 
solutions allowed by the symmetries and dynamics of the system must be expressed in 
the wavelet basis, i.e. they must be a complete set. Kaiser 50 has formulated relativistic 
electrodynamics in terms of a wavelet based on analyicity and coherent states. D’Aranio 
and DeFacio 49 have formulated a class of squeezed states of quantum optics and provided 
a general inversion structure for a general density operator in terms of a suitable window, 
which was expressed in terms of a general basis of observables. Han, Kim and Noz 51 have 
presented a study of the relativistic phase space of light using compactly supported wavelet 
windows to clarify the relation between photons and light waves. Since there are many 
wavelet windows for any given problem, it is important to know which ones are natural 
and, if possible, if a “best one” exists. 

In this paper, the application of wavelets and state squeezing, Bogoliubov transfor- 
mations to the two statistical mechanics problems: E (q, u?) for water in Section 2 and 
bubble formation will be described in Section 3, and our conclusions will be presented in 

Section 4. 
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2. STATISTICAL MECHANICS OF WATER, THE DIELECTRIC RESPONSE 

In collaboration with Professor G. Vignale, we have published studies of a classical 
water-like polar fluid whose properties are chosen to mimic water. The starting point is the 
number- conserving relaxation- time approximation to the Boltzmann transport equation. 
Up to the first order deviation from equilibrium it has the form 



where overdots are time rate of changes, q a — p a — — fo = e ^ l Z is the equilibri- 
um Boltzmann distribution function. T is the kinetic energy of the molecule and includes 
translations, rotations and small amplitude vibrations, and Z is the partition function. 
f\ is the first order correction term to /o by a perturbed applied potential. The static 
response of the electric susceptibility x(w = 0) is a 3 x 3 real matrix-valued quantity given 
by -Pfo(q a )- . 

Upon determining the dynamical self-response function self-consistently from /j 


X.(“>) = (1 + i/wr)x 0 (w + i / t ) 


1 + — xV + i/r) {x°(u> = 0)} 
t or 1 } J 


(ii) 


where the Kubo fluctuation dissipation theorem for the single particle self-electric suscep- 
tibility 


X°M = x(0) [1 + *«G(«)] 


( 12 ) 


relates the electric susceptability to the van Hove correlation function G(u). The total 
response from the Kubo theorem is 

x(w) = x.(w) [i - *(9)x.H] _1 » ( 1 3 ) 

where 'k is the local field factor which describes the coherent many particle interactions. It 
is the analog of Lorentz’s local field factor, in modern many-body theory. We have explored 
several choices and found Wei and Patey 17 to provide the best agreement with experiment. 
Indeed, various computer intensive molecular dynamics and Monte Carlo calculations gave 
er(0,0) as 25-71 instead of the experimental value of 90.8 ± 3.2! The large experimental 
uncertainty is dominated by the spontaneous dissassociation H^O ?=* H + + OH~ , which 
at 20° C yields an ion concentration of 8.33 x 10" 8 moles /liter. The high mobility of the 
H + species allows it to form a positively charged surface layer at the negative electrode 
which can distort the measured value of the electric permittivity, e(u;). By taking a limit 
as the moments of inertia approach zero, the Stockmayer fluid calculations were found 
comparable to the best computer studies. The causal dielectric permittivity values are 


er(9» [1 + 47rxr(<7>)] 

1 

CL(q,u) 


and 


= [1 - 4ffXi(£ “)] 


(14b) 



where (T, L) denote the transverse and longitudinal components of the dielectric tensor e 
and the electric susceptability tensor *. The directions (T, L) are taken from the direction 
of the applied (vacuum) electric field E. Remark: A wave propagation or scattering 
experiment is required to measure the q dependence of n and a, in contrast to the ac 
capacitance bridge which measures only the angular frequency, «, dependence. A graph 
from ref. (18) of the real part of the index of refraction will be shown later. 

There is a collective mode in e L which Lobo, Robinson and Rodrigues 9 first suggested 
and Pollock and Alder 10 found in their simulation of a Stockmayer fluid model of water. 
It occurs in our model both in the symmetric rotor “waterlike” case 18 and in the realistic 
asymmetric rotor case. 19 

Considering the water molecule as a two-level system for simplicity, the composite Boson 
creation and destruction operators which create and destroy water molecules ( recall that 
this analysis is at frequencies that are far below the breakup threshold energies) are written 
as (««, «£) where their directions are those of the electric dipole moment, £, of the molecule 
and <*,/? = 1,2 label which level. The collective dipolaron boson is created and destroyed 
by the Bogoliubov transformed operators (, b a , 6J). The frequency of the dipolaron is 
complex- valued 


<^>a Qqa ^T^a (lb) 

and lies in the lower half w-plane because the response function l/e L (q,u) is causal. The 
frequency in eq. (15) is determined from 


Re [e L (q, ft, Q )] = 0 


and the width from 


(16) 


p _ |/m[e L (},n,)]| 

’“ = I ML 


(17) 


provided that Cl qa » r ?a . The peak and its width are shown in Fig. 1, where one c an 
see that Q, qa >> r ?a is satisfied. The mechanism of this collective mode is high frequency 
coherent oscillations of the electric dipoles in self-consistent fields in the liquid. In order 
for the mode to be longitudinal, they oscillate 180° out of phase. Both Q qa and T qa were 
shown by Kim, et al. 1& , to be almost dispersion-free (they only change by 1% over the 
entire range of validity). This mode depends on q = \q\ being small but non-zero so it 
cannot be measured from any ac capacitance experiment such as those in ref. (3) and must 
await a light or neutron scattering experiment. Such an experiment would verify or falsify 
the model formulated in refs. (18,19). 
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(■> (b) 
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Fig. 1. The dipolaron collective mode in water. 

Next, the reader is reminded that eq. (10) includes translations, rotations and vibrations 
in the p a (§jr) tenns - It is well known 52,53 that the renormalization group is based upon a 
scaling invariance. By expanding the single particle distribution function /i ( •) in a wavelet 
series 

/.(•) = £<VM-) (is) 

ji 

where (•) can stand for either the translation x, rotation (9, <f>) or vibrations ( 91 , 92 , 93 ) i n 
the real index of refraction n(u>). The three normal modes for vibration are given by group 
theory and are used in our calculation of n(u). In Fig 2 the frequencies up to 10 10 Hz 
are dominated by the translations, from 10 10 ~ 10 n Hz at (A, A') the collisions dominated 
Debye relaxation time dominates and at 10 l3 Hz the free rotor peak occurs at (Bi, £'). 
The free rotor peak at ( B ') in our theory is much too large, so work is underway to try 
to bring this into better agreement with the experiment. Between 10 13 and 10 u Hz the 
collective mode (C') occurs with coherent oscillations of the dipoles as its mechanism. The 
small optical peaks at (D U D 2 ) fail to resolve the two higher frequency modes in (D 2 ) but 
we think that future experiments will. Thus, the wavelet phase space allows us to analyze 
the Boltzmann transport equation with a non-trival interaction. 
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Fig. 2. The real index of refraction of water. 


3. BUBBLE FORMATION IN WATER BEING INSONIFIED BY ULTRA- 
SOUND 

The idea behind this calculation is to adapt the idea of Kaup to a solitary wave model 

by generalizing the work of Hammer, Shrauner and DeFacio where <j> is a 1 © 1 fluid density 
field 


L 


~2 M *) 2 + (^<) 2 



(19) 


with u$ = A/B where A, B are real and positive coupling constants and A is the strength 
of the negative mass, B is the strength of the positive nonlinear quartic interaction to three 
space dimensions with a radial time independent potential. In 3 ® 1 space-time dimensions 


L f = ~ [( v </>) +^< 2 

La = -\[(Vp) 2 -p[ 
Li = gtfp 2 + p<f>] 



(20a) 

(20b) 

(20c) 


where L F is the Lagrangian of the fluid, L A is the acoustic Lagrangian for the sound wave 
and Li is the interaction Lagrangian. If p. ~ 0 and the sound waves are small enough in 
amplitude to be linear (see eq. (30a)) the field equations become 


- (A <f> -(- <f > tt ) - A<f> + B(f> z + gp 2 <j> = 0 , (21a) 

- (Ap-ftt) - Aip + Bxp* + g<f> 2 p = 0 (21b) 

The key role played by the nonlinear “potential” which acts as a convection term is our 
justification for calling this approach a strong convection model. The mechanism is that 
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the constant sound wave in the fluid produces gas under pressure as a new phase in the 
fluid. This is a non-equilibrium, open system and the gas density can be described by one 
Bogoliubov transformation and the non-zero temperature, 0 °C to 100°C, can be obtained 
from another Bogoliubov transformation. The entropy of this situation was discussed by 
DF, VN and Professor Brander 63 of the Institute for Theoretical Physics at Chalmers in 
the Gothenberg. The (. = 0, radial solutions to eqs. (21a, b) are Jacobi elliptic functions 
with modulus k u 0 < h < 1, <5 0 i is a density defect where the local mass density is a 
minimum so that it is a site for cavitation and (u,uqi) are the amplitudes for the density 


and pressure wave, respectively. 

Next a time-independent mean field theory for the density <f> 
about the Jacobi elliptic functions according to 

is formed by linearizing 


gp 2( f> -* 9 u h [Sn(Ki (r + <5 0 i) \h)] 2 < f > , 

(22) 

and 

(23) 

B<t> z —> B < <f> 2 > <j> -> Bu 2 [Sn (k (r + 8 ) |fc)] 


This should be valid, provided it is not applied too close to a caviton ‘explosion’, and leads 
to the linearized field eqaution for the density <f> 

A <f> + <t>u + Vmft('>') ( I> = 0 ^ (24) 


where 

V MFT (r) = A- Bv? [Sn (/c (r + 8) |fc)] 2 + gu\ x [Sn (k x (r + <5 0 i) |*i)]* 


Separation of variables using 


and 


<f> = R(r)Y tm (e^)e-^ 


JJ(r) = 


u(r) 


(25) 

(25a) 

(25b) 


reduce eq. (25) to 


Urr + 


Ymft ~ 






u = 0 


(26) 


This ODE satisfies Dirchlet zero boundary conditions on the dense domain H ( ) where 
R x + is the positive half-line and ff 2 is the Sobolev space of L 2 functions whose first two 
derivatives exist and are X 2 . Thus, eq. (26) can be solved numerically using the Numerov 
algorithm. Then the density is constructed from the relation 

*) = £ c q i m R q t(r )Yt m {9, 4>)e . (27) 


From eq. (27) the sound wave has modulus k t ~ 0. One thing about this strong convection 
model which seems new is the inclusion of £ > 1 multipoles which give deviations from 
spherical symmetry. The calculus of variation “proof” that bubbles are spherically shaped 
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only holds for thermal and fluid equilibrium, whereas the continous quasi-monochromatic 
sound wave makes this an open, non-equilibrium system. In Fig. 3 Vmft is plotted and 



Fig. 3. An example of a mean field theory potential, Vmpt- 



^*8* 4. The water density in a 3 mm cell calculated from the mean field theory 

potential in the previous figure. 

in Fig. 4 a bubble is shown which was found in the density, using the convection potential 
of Fig. 3. The parameters are consistent with those of water with a sound wave amplitude 
which is physically realizable. The bubble is centered at 6 = 0.75mm; and has radius 
fQ 300/x (/x 10 m). The detection of bubble interfaces using wavelets is under study. 

Since some authors have questioned the safety of ultrasound in tissue, 31,34 it is important 
to understand bubble formation. In addition, the Mie scattering from bubbles adds spatial 
dispersion to the effective index of refraction. 
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4. CONCLUSIONS 

Two different problems in non-equilibrium statistical mechanics were discussed. 

(i) the dielectric response of water, e(<f, u>), and 

(ii) bubble formation in water under insonnification using Bogoliubov transformations: 
state squeezing maps and wavelets. 

The multi-resolution structure of phase space permits the sort of analysis which the 
renormalization group has provided for lsing models and lattice gas models and little else. 
Thus, “squeezed states” are useful for more than quantum optics. 
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